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II. “ On a Meteoric Stone found at Makariwa, near Invercar¬ 
gill, New Zealand." By G. H. F. Ulrich, F.G.S., Pro¬ 
fessor of Mining and Mineralogy in the University of 
Dunedin, N.Z. Communicated by Professor J, W. Judd, 
F.R.S. Received December 14, 1892. 

(Abstract.) 

The specimen described in this memoir was found in the year 1879 
in a bed of clay which was cut through in making a railway at In¬ 
vercargill, near the southern end of the Middle Island of New 
Zealand. Originally this meteorite appears to have been about the 
size of a man’s fist, and to have weighed 4 or 5 lbs., but it was 
broken up, and only a few small fragments have been preserved. The 
stone evidently consisted originally of an intimate admixture of 
metallic matter (nickel iron) and of stony material, but much of the 
metallic portion has undergone oxidation. Microscopic examination 
of thin sections shows that the stony portion, which is beautifully 
chondritic in structure, contains olivine, enstatite, a glass, and prob¬ 
ably also magnetite, and through these stony materials the nickel 
iron and troilite are distributed. The specific gravity of portions of 
the stone was found to vary between 8*31 and 3*54, owing to the un¬ 
equal distribution of the metallic particles. A partial chemical 
examination of this meteorite was made by the author and Mr. James 
Allen, but the complete analysis has been undertaken by Mr. L. 
Fletcher, F.R.S., of the British Museum. The analysis, which when 
finished will be communicated to this Society, has gone so far as to 
show that the percentage mineral composition of the Makariwa 
meteorite may be expressed approximately by the following numbers : 
nickel iron 1, oxides of nickel and iron 10, troilite 6, enstatite 39, 
olivine 44 


III. “ On Operators in Physical Mathematics. Part I.” By 
Oliver Heaviside, F.R.S. Received December 15, 1892. 

Connexion between a Flux and a Force through an Operator. 

1. In the investigation of physical questions we often have ta 
answer such a question as this : Given a force /, a function of the 
time, acting at one place in a connected system, find the effect F, of 
some given type, produced by the force at its own or some other 
place. Or it may be that it is not an impressed force that is given, 
but displacement of some kind. Or, in order to produce mathe- 
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matical simplicity, we may liave a space-distribution of force or of 
displacement given, whose effect is required. 

To answer the question, we may investigate the general differential 
equation of the system, find its solution (series, integrals, &c.), and 
then introduce special values of constants or of functions to limit the 
generality of the problem, and bring the solution to satisfy the re¬ 
quired conditions. Details may differ according to circumstances, 
but this may serve to describe the usual process. 

2. There is, liowever, a somewhat different way of regarding the 
question. We may say that we have no special concern with, the 
general solution which would express the disturbance anywhere due 
to initial energy throughout the system; but that we have simply a 
connected system, a given point (for example) of which is subjected 
to impressed force, communicating energy to the system, and we 
only want to know the effects due to this force itself. Since, 
therefore, the connexions are definite, we must have some definite 
connexion between the “flux ” P and the “ force ”/, say 

P - Yy, (1) 

where Y is a differentiating operator of some kind, a function of djdt , 
the time-differentiator, for instance, when the connexions are of a 
linear nature. Here / is some given function of the time, and Y 
indicates the performance upon / of certain operations, whose result 
should be to produce the required function P. 

3. An important point to be noted here is that there is, or should 
be, no indefiniteness about the above equation. The operator Y 
should be so determined as to fully eliminate all indeterminateness, 
and so that the equation contains in itself the full expression of the 
connexion between the force and the flux, without any auxiliary con¬ 
ditions, or subsequent limitations, except what may be implicitly 
involved in the equation itself. 

Determinateness of a Solution through the Operator. 

4. But as soon as we come to distinctly recognize this determin¬ 
ateness of connexion, another point of important significance presents 
itself. It should be possible to find P completely from / through the 
operator Y without ambiguity and without external assistance. 
That is to say, an equation of the form (1) not only expresses a 
problem, but also its solution. It may, indeed, not be immediately 
interpretable, but require conversion to some other form before its 
resultant meaning can be seen. But it is, for all that, a particular 
form of the solution, usually a condensed form, though sometimes it 
may be of far greater complexity than the full ordinary solution. In 
this respect the nature of the function / is of controlling importance. 
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We need not assert that the determinateness of E from equation (1) 
is true for all forms of the function Y that may be written down 
arbitrarily ; but that it is true in the forms presenting themselves in 
dynamical problems seems to be necessitated. 

5. We have, therefore, presented to us the problem of solving this 
equation for any particular form of Y that occurs. This may be 
very easy and obvious, or it may be excessively difficult and obscure. 
In the latter case it may be so merely because we do not know how 
to do it. Then we should find out. As our argument is that Y finds 
E from/definitely, there should be definite rules for the manipula¬ 
tion of the operator Y, or of the expression Y/, for its conversion to 
the form of an ordinary mathematical function, which will be the 
solution in the usual sense, freed from differentiating operations. We 
may find how to work by experiment. For, if two different methods 
lead to different results, one of which we find to be correct by in¬ 
dependent tests, we can safely assert that one of the methods was 
partly wrong, whilst the other may have been wholly correct. So by 
practice we may come to know something about it. 

6. Again, the function Y, regarded as an algebraical function, may 
admit of different forms of expression. These are algebraically 
equivalent, but to what extent they may be equivalent in their ana¬ 
lytical aspects—for instance, one series involving differentiations 
equivalent to another involving integrations, and leading to results 
which are either identical or equivalent—cannot be safely said 
beforehand. It is, in its generality, a rather difficult and obscure 
matter. In special cases I find that forms of Y which are algebrai¬ 
cally equivalent are also analytically equivalent; but I have not 
succeeded in determining the amount of latitude that is permissible 
in the purely algebraical treatment of operators. No doubt there are 
definite limitations, but they have to be found. I have, however, 
extensively employed the algebraical treatment experimentally,* sub¬ 
ject to independent tests for guidance. It proved itself to be a 
powerful (if somewhat uncertain) kind of mathematical machinery. 
We may, for example, do in a line or two, work whose verification 
by ordinary methods may be very lengthy. On the other hand, the 
very reverse may be the case. I have, however, convinced myself 
that the subject is one that deserves to be thoroughly examined and 
elaborated by mathematicians, so that the method may be brought 
into general use in mathematical physics, not to supplant ordinary 
methods, but to supplement them ; in short, to be used when it is 
found to be useful. As regards the theory of the subject, it is in¬ 
teresting in an unusual degree, and the interest is heightened by the 
mystery that envelops certain parts of it. 

# The reader will find examples in my ‘ Electrical Papers/ vol. 2, of the treat¬ 
ment of irrational as well as rational operators. 
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Electromagnetic Operators. 

7 . Perhaps the best way of beginning the subject, to obtain a 
good idea of the nature of the operators and the advantages 
attending their use, is through the theory of a connected system 
of linear electrical conductors. For the electrical equations seem 
to be peculiarly fitted for the illustration of abstract dynamical 
properties in a clear manner, even when quite practical electro¬ 
magnetic arrangements are concerned. We know that we may, by 
the application of Ohm’s law to every conductor (or to circuits of 
conductors), express the steady current G n in a conductor n due to an 
impressed force e m in a conductor m by an equation 


0 % — Y mnCmt (2) 

where Y mn is some algebraical function of the resistances of the con¬ 
ductors—usually of all the resistances, although in special cases it 
may become independent of the values of some of them. Now, 
suppose it is not the steady current that is wanted, but the variable 
current when e m varies. The answer is obviously given by the same 
equation when the function Y involves only resistances; that is, when 
there is no storage of electric or magnetic energy, so that Y is a con¬ 
stant not involving djdf. Then the flux and the force keep pace 
together, and their ratio does not vary. It is, however, less obvious 
that the same equation should persist, in a generalized form, when 
every branch of the system is made to be any electromagnetic 
arrangement we please which would, in the absence of its connexions 
with the rest of the system, be a self-contained arrangement. To 
obtain the generalized form of Y mn we have merely to substitute for 
the resistances concerned the equivalent resistance operators. That 
is, instead of Y = HO, where Y is voltage, G current, and R resist¬ 
ance, we have an equation Y = ZO in general for every conductor, 
where Z is the resistance operator appropriate to the nature of the 
conductor, which may be readily constructed from the electrical 
particulars. These Z’s substituted in Y nm in place of the R/s make 
equation (2) fully express the new connexion between the flux 
and the force e m . There is much advantage in working with resis¬ 
tance operators because they combine and are manipulated like 
simple resistances. Of course (2) is really a differential equation, 
though not in the form usually given. To make it an ordinary 
differential equation we should clear of fractions, by performing such 
operations upon both sides of (2) as shall remove denominators and 
all inverse operations. It is then spread out horizontally to a great 
length (usually) and becomes very unmanageable. Also, we lose sight 
of the essential structure of the operator Y. 

8. By arrangements of coils and condensers in our linear system 
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we may construct an infinite variety of resistance operators, and of 
conductance operators, such as Y mn above. They are, however, always 
algebraical functions of p, and are finite. If expanded, equation (2) 
always becomes an ordinary linear equation of a finite number of 
terms. But if we allow conduction in masses, or dielectric displace- 
ment in masses (with allowance for propagation in time), the finite 
series we were previously concerned with become infinite series. 
This, at first, appears a complication, but it may be quite the reverse, 
for an infinite series following an easily recognized law may be more 
manageable than a finite series. Still, however, the equivalence to 
ordinary differential equations persists, provided our arrangement is 
bounded. But when we remove this restriction, and permit free dis¬ 
sipation of energy in space (or equivalently), another kind of operators 
comes into view. The complexity of the previous, due to the reaction 
of the boundaries, is removed; simpler forms of operators result, and 
they do not necessarily admit of the equations taking the form of 
ordinary differential equations, as they may be of an irrational nature. 
This brings us necessarily to the study of generalized differentiation, 
concerning which, more presently. 


Operators admitting of Easy Treatment . 

9. In the meantime, notice briefly some of the ideas and devices 
that occur generally in the treatment of operators. First of all, we 
may obtain the steady state of F due to steady /, when there can be 
a steady state of F, by simply putting p = 0 in the operator Y con¬ 
necting them, p meaning djdt. Even when there is no resultant 
steady state of the flux, as when reflections from a boundary continue 
for ever, the term F = Y 0 /has its proper place and significance. 

Next, we may notice that if the form of Y should involve nothing 
more than separate differentiations, as in 

F = (a + bp + cp*+....)f, (3) 

then all we have to do is to execute the differentiations to obtain F 
from /. When / is a continuous function, this presents nothing 
special. When discontinuous, however, a special treatment may be 
needed. 

In a similar manner, there may be only separate integrations or 
inverse differentiations indicated in Y, as when 

F = (a - f bp- 1 + ep~*+ . ...)/. (4) 

Since / is a definite function of the time, so are its successive time- 
integrals. In this case, / may be discontinuous, and yet present no 
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difficulty. Suppose it is zero before and constant after the moment 
t — 0. Then we shall have 

in 

V= (5) 

P. 

A combination of direct and inverse operations, which frequently 
occurs in the theory of waves, is exemplified in 

F = c~P r l v (a -f bp~ l -f cp~ 2 + (0) 

Here we may perform the integrations first, getting the result 
fp(f) say, and then let the exponential operate, giving, by Taylor’s 
theorem, 

F = 0(*-r/v). (7) 

Or we may let the exponential operator work first, and then perform 
the integrations. This may be less easy to manage, on account of 
the changed limits. 

Two important fundamental cases, which constitute working 
formulae, are 

and F =-£-/, (8) 

p — a p+a 

with unit operand, that is, /=0 before and constant after t — 0. 
Here we may expand in inverse powers of y>, getting, in the first 
case 

F = (l + ap~ l -f- a 2 p~ 2 + •. .. ) = e®*, (9) 

and in the second case e”®*. The latter expresses the effect of a 
unit impulse in a system having one degree of freedom, with friction, 
as when an impulsive voltage acts upon a coil. 


Solutions for Simple Harmonic , Impulsive , and Continued Forces. 

10 . A very important case, admitting of simple treatment, occurs 
when the force is simple periodic, or a sinusoidal function of the 
time. It may happen that the resulting state of F is also sinus¬ 
oidal. For this to occur, there must be dissipation of energy, to 
allow the initial departure from the simple periodic state to subside. 
We then have_p 2 = — n 2 applied to F as well as f, where n)2ir is the 
frequency ; so that the substitution of ni for p in Y brings equation 
(1) to the form 

F = (Y 0 +Y iij/ = (Y 0 + Y in -'p)f, (10) 

where Y 0 and Yi are functions of n 2 . We now find F by a simple 
direct operation. This case is so important because its application 
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is so general, and its execution usually presents no difficulties, whilst 
the interpretation of the result may be valuable and instructive 
physically. 

A continued constant force of unit strength, commencing when 
t — 0, may be represented by 


F = 



sin nt 
n 



(ii) 


using a well-known integral. We may apply this to equation (1), if 
desired, and obtain a particular form of solution. And from (11) we 
see that a unit impulse is represented by 


P 


WJ 0 


cos nt . dn 


( 12 ) 


acting at the moment t = 0. This is, of course, the basis of Fourier’s 
theorem. But, instead of the application of the fully developed 
Fourier’s theorem, it is more convenient to use (12) itself. Thus, 
when / 0 is an impulse acting' when t = 0, we have the equation 

F = Yp/o, 

p/o representing the force. So, by (12), 

F=Af°Ycos nt.dn (13) 

if J 0 

gives us a particular form of the solution arising from an impulse. 
Take p = ni in Y to convert the quantity to be integrated to an 
algebraical form. 

Since a continuously varying force may be represented by a suc¬ 
cession of infinitesimal impulses, we see that a single time-integration 
applied to (13), f 0 being then a function of the time, gives us a form 
of solution of the equation F = Y f, for any kind of / and Y that can 
occur. It is, however, a theoretical rather than a practical form of 
solution. For it usually happens that the definite integral is quite 
unamenable to evaluation. The same may be often said of the solu¬ 
tion (13) for an impulse, and in such cases it may be questioned 
whether the form F = Y/itself is not just as plain and intelligible. 
In fact, in certain cases, a very good way to solve (or evaluate) a 
solution in the form of a definite integral is to undo it, or convert it 
to the symbolical form F = Y f, and then solve it by any way that 
may be feasible. Nevertheless, it is interesting to know that we may 
have a full solution, and the definite integrals are sometimes practic¬ 
ally workable, or may be transformed to easier kinds. 
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11 . There is also the method of partial fractions. It is not always 
applicable, and is especially inapplicable when the removal of 
bonndaries drives the roots of the determinantal equation into con¬ 
tiguity. But the application is very wide, nevertheless. Put 
Z = Y” 1 ; then the solution of F = Y/, when / is constant, starting 
when t — 0 , is 

ai) 

where Z 0 is the steady Z, got by taking p = 0 in Z, and the summa¬ 
tion ranges over the roots of the equation Z = 0, considered as an 
algebraical equation in p . That is, p is entirely algebraical in (14). 
Similarly, the effect of an impulse f 0 is represented by 




ePt 

dZijdp ’ 


(15) 


and from this again, by time* integration, we can obtain an expression 
for the effect due to any varying /, which may be quite as un¬ 
manageable as the previous definite integral for the same. On the 
other hand, (14) and (15) furnish the most direct and practical way 
of investigating certain kinds of problems, whether there be but a 
few or an infinite number of degrees of freedom. This method is 
the real foundation of all formulae for the expansion of arbitrary 
functions in series of normal functions. For, find the impressed 
forcive that would keep up the arbitrary state. We may then apply 
the above to every element of the forcive to find its effect, and by 
integration throughout the system get the arbitrary functions ex¬ 
panded in normal functions. Or, without reference to impressed 
force, find the differential equation connecting any element of the 
initial state and the effect it produces later. It will be of a form 
similar to our F = Y/, and it may be similarly solved by a series, 
which contains the expression of the expansion of the initial state in 
the proper functions. 

Or we may investigate the normal functions themselves, and em¬ 
ploy their proper conjugate property to obtain the expansion repre¬ 
senting any initial state. But this method does not apply very 
naturally to equations of the form we are considering. 


Decomposition of an Operator into a Series of Wave Operators. 

12 . There is also another method which contrasts remarkably with 
the previous, viz., to decompose the operator Y into a series of other 
operators of a certain type expressing the propagation of waves. This 
is best illustrated by an example. Suppose the question is, given an 
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impressed force acting at one part of a long telegraph circuit, find the 
effect produced. One way would be to first find the effect due to a 
simple periodic force ; from this the effect of an impulse follows ; and 
from the latter the effect due to any/. A second way is by means of 
the normal functions, either through the conjugate property or by 
partial fractions. Lastly, we may decompose the operator Y into 
operators of the form which would exist were the circuit infinitely 
long, so that the effect of terminal reflections and absorptions does not 
appear. Say we have 

F = (Y 0 + Y 1 + Y, + ....)/. (16) 

Then F = Y 0 f will represent the initial wave from the source /, 
whilst the rest will express the succeeding reflected waves from the 
terminations of the circuit. The operators Y 0 , &c., may be all of the 
same type, so that it suffices to solve F — Y 0 /, that is, convert it to 
an ordinary algebraic functional form, to obtain that form of the 
complete solution which has the greatest physical meaning, inasmuch 
as it shows in detail the whole march of F in terms of /. So does the 
solution in terms of normal functions, but not immediately, because 
the successive waves are expressed in the form of an infinite series of 
vibrating systems. Their resultant effect cannot be seen at once. 
We might, indeed, almost say that the form of solution in successive 
(or simultaneous) waves was the solution, being of the most explicit 
nature. Should, however, the impressed force be of a distributed 
nature, of the type suggested by a normal function, for example, 
then clearly it is the expression in terms of waves that becomes com¬ 
plex and unnatural. We also see that, although a direct transforma¬ 
tion from one form of solution to another may be wholly impracticable 
algebraically, yet it may be readily carried out through the function 
Y as intermediary. 


Treatment of an Irrational Operator. Solutions in Ascending Series. 


13 . The above general remarks are necessarily very sketchy. Some 
of the matters mentioned may be returned to, but the object of the 
preceding is merely to prepare the mind of the reader for the more 
transcendental matter to follow. Let us now consider how to treat 
irrational operators directly, without the assistance of definite integ¬ 
rals. The first form that presented itself to me was that exhibited 

fey 


/K + S/Ai 
\R + Lp) ’ 


( 17 ) 


where p is d/dv and R, S, K, L are constants. It occurs in the theory 
of a submarine cable or other telegraph circuit, and in other problems. 
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R and L are themselves differentiating operators of complicated form 
in general, or, more strictly, R + Lp is a resistance operator, say R", 
of very complex form. Bat it is quite sufficient to take the form 
R + Lp, where R is the effective resistance and L the inductance per 
unit length of circuit. S and K mean the permittance and leakage 
conductance per unit length. 

Now we may readily obtain the simple periodic solution out of (17), 
by the before-mentioned substitution p = ni; and in doing so we may 
use the general operator R", for that will then assume the form 
R + Lp. From this solution a wholly uninterpretable definite integral 
can be derived to express the effect of an impulse or of a steady im¬ 
pressed force. The question was, how to obtain a plain understand¬ 
able solution from (17) itself to show the effect of a steady force. 
To illustrate, we may here take merely the case in which K = 0, whilst 
R and L are constants, because the inclusion of K (to be done later) 
considerably complicates the results. We have then to solve 


F = 



(18) 


where a is a constant and p — djdt . The operand is understood to 
be unity, that is, / — 0 before and =1 after t = 0. It is needless to 
write unit operands, and it facilitates the working to omit them. 
Now, the first obvious suggestion is to employ the binomial theorem 
to expand the operator. This may be done either in rising or in 
descending powers of p. Try first descending powers, since by ex¬ 
perience with rational operators we know that that way works. We 
have 


F = (l + ap- 1 )"* 


j a l. 3 a 2 

* p 2 2 |2 p 2 


1.3.5 a 3 


213 


P 


+. (19) 


The integrations, being separated from one another, can be imme¬ 
diately carried out through p~ n = t n j\n, giving the result 


1.3 / at \ 2 1.3.5/ at\ 


or, which is the same, 



( 20 ) 


( 21 ) 


where I 0 is the well-known cylinder function, Now, that this re¬ 
sult is correct may be tested independently, viz., by its correctly 
satisfying the differential equation concerned and the imposed con¬ 
ditions. We therefore obtain some confidence in the validity of the 
process employed. 
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In a precisely similar manner we may show that 
P 


P— a . 


2 T I at 

= e 2 Io o 


( 22 ) 


Further modifications are confirmatory. Thus, by making use of 

/ (p) c at — e*‘f(p + a ), (23) 

we can shift e at and similar functions back and forth. Also using the 
stock formulas 


P — e at P — e ~at 

p-a~ ’ p+a~ ’ 

we have the following transformations, 



(23a) 


(24) 


In this we may use the result (22), and so come round to (21) again. 
From the above we see that 


I 0 («0 — 


e at 


P 

2a+p i 


= (T at 



(25) 


and further, by shifting the exponentials to the right, to make them 
the operands (instead of t °), 


I 0 (ai) = 


p~ a 
\P + a 


(p±a\* 

xP~~ a ) 


e at ; 


(26) 


and now further again, by employing (23a) in place of the exponen¬ 
tials, we obtain 


Io (at) = 


r p-~a\ k p _ p 

JP + ®/ p—a ~~ (p 1 — a 3 )^’ 


(27) 


which is an entirely different kind of operator, since the square of p 
occurs under the radical sign, instead of the first power. But (27) 
may be readily tested and found to be not wanting. For expand by 
the binomial theorem, thus, 


_ P _ 

(p 2 ~a 2 y 



i+i 


a 3 1.3 a 4 

p 2 2 3 |2 _p 4 


(28) 


This may be immediately integrated, giving as result the series 

T (n A _ , , (« 0 3 . (aty , 

\ at ) — 1 + — - ’ 


(29) 
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the well-known formula for I 0 (a£) in rising powers of the square of 
the variable, as required. 


Transformation to a Descending Series . 

14 . There is such a perfect harmony in all the above transforma¬ 
tions, without a single hitch, that you are tempted at first to think that 
you may do whatever you like with the operators in the way of alge¬ 
braical transformation. There is a considerable amount of truth in 
this, but it is not wholly true. I shall show later some far more 
comprehensive and surprising transformations effected by simple 
means. At the same time I should emphasize the necessity of 
caution and of frequent verification, for no matter how sweetly the 
algebraical treatment of operators may work sometimes, it is subject 
at other times (owing to our ignorance) to the most flagrant failures. 

But in the above we only utilized one way of effecting the binomial 
expansion. There is a second way, viz., in ascending powers of the 
differentiator. The two forms are algebraically equivalent so far as 
the convergency allows, but we have, so far, no reason to suppose 
that they are analytically equivalent. But on examination we find 
that they are. Thus, using the first of (25) and expanding, we 
get 


Io (at) 


c at 


p v 

2a+p) 




= ^(i 

1 4a 2 



(30) 


Here the operand is t° or unity. Or we may make it ( pj2a )* if we 
please. If we know its value, as a function of t , the rest of the work 
is easy, as it consists merely of differentiations. But nothing that 
has gone before gives any information as to the meaning of p*, let 
alone its value. We may, however, find it indirectly. We may 
prove independently that when at is very big, 1 0 (at) tends to be re¬ 
presented by e at (2 7r at)~K From this we conclude that the value of 
p 4 must be (7 rt)~K Then (30) becomes 


I„(ai) = 6 W 1-A+— (- 
v J ' 4 a 2 \4a 


(30a) 


1(2 7 rat) k ’ 

and now performing the rest of the differentiations, we arrive at 

Io(o/) = (i^so* { 1+ i + |i(8^r + i^ 3 5p + ••••}' (31) 


which, on test, is found to be equivalent to the ascending series (29). 
Of course only the convergent part of the series can be utilized for 
VOL. lil. 2 M 
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calculating the value of the function. It is, however, the series for 
practical use when at is big enough to make calculation by the con¬ 
vergent ascending series very lengthy. Stop when the convergency 
of (31) ceases. The result will be too big. Leave out the last 
counted term, and the result is too small. Counting only half the 
last convergent term, the result is nearly right, being a little too 
big. There seems no possible way of hitting the exact value. But 
still, when at is big, we can get quite close enough, to four or 
more figures, or any other number we please when at is sufficiently 
increased. 


Fractional Differentiation. 

15 . Knowing in the above manner _pt, the values of pi, pi, &c., follow 
by complete differentiations. But although, on the basis of the above, 
a considerable amount of work may be done, and extensions made, 
yet it is desirable to stop for a moment. For the whole question of 
generalized differentiation is raised. The operator pi presents itself 
in analogous problems, along with pi, &c. We want a general method 
of treating p n , when n is not confined to be integral. Notice, however, 
in passing a remarkable peculiarity of the above investigation. If 
we had put L = 0 in (17), as well as K = 0, we should have had the 
form Y — pi to consider at the beginning, with no evident means of 
treating it. By taking, on the other hand, a more general case, as we 
did, we avoided the fractional differentiation altogether, and easily 
obtained a convergent solution, viz., (21), through (18), (19), and 
(20). It is not always that we simplify by generalizing. 

The sum total of the whole information contained in my mathe¬ 
matical library on the subject of generalized differentiation is con¬ 
tained in the remark made on p. 197 of the second part of Thomson 
and Tait’s ‘Natural Philosophy/ paragraph ( n ), relating to the 
process by which spherical harmonics of any degree may be derived 
from the reciprocal of a distance :—“ The investigation of this gene¬ 
ralized differentiation presents difficulties which are confined to the 
evaluation of P 5 , and which have formed the subject of interesting 
mathematical investigations by Liouville, Gregory, Kelland, and 
others.” 

I was somewhat struck with this remark when I first read it, in 
trying to plough my way through the fertile though rather heavy field 
of Thomson and Tait, but as the subject was no sooner mentioned than 
it was dropped, it passed out of mind. Nor did the absence of any 
reference to the subject in other mathematical works, and in papers 
concerning mathematical physics generally, tend to preserve my re¬ 
collection of the remark. Only when the subject was forced upon my 
attention in the above manner did I begin to investigate it, and not 
having access to the authorities quoted, I was compelled to work it out 
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myself. I cannot say that my results are quite the same, though there 
must, I think, be a general likeness. I can, however, say that it is a very 
interesting subject, and deserves to be treated in works on the Integral 
Calculus, not merely as a matter concerning differentiation, but 
because it casts light upon mathematical theory generally, even upon 
the elements thereof. And as regards the following brief sketch, 
however imperfect it may be, it has at least the recommendation of 
having been worked out in a mind uncontaminated by the prejudices 
engendered by prior knowledge acquired at second hand,. I do not 
say it is the better for that, however. 


Differentiation Generalized „ 

16. The question is, what is the meaning of if v signify d\dx , 
when n has any value ? This is, no doubt, partly a matter of con 
vention; but apart from all conventions, there must be fundamental 
laws involved. How observe that the effect of a whole differentiation 
V upon the function x' 1 is to lower the degree by unity. This applies 
universally when n is not integral. When it is integral, there seem 
to be exceptions. But we can scarcely suppose that there is a real 
breach of continuity in the property. We also observe that a whole 
differentiation v multiplies by the index, making y yx n = nx n ~ l ; and 
again there are apparent exceptions. How the first thing to do is to 
get rid of the exceptions. Hext, the obvious conclusion from one v 
lowering the index by unity, v 2 by two, and so on, is that y >l lowers 
the degree n times, whether n be integral or fractional. Further, since 

x n X n-1 

v r — i— i 9 

I n \n — 1 


when n is positively integral, and | n is the factorial function 
1 .2.3 .. .. w; and, similarly, 



(33) 


whatever positive integer n may be, it is in agreement with the 
previous to define generalized differentiation by the last equation, for 
all values of n , provided we simultaneously define |rp to be given by 


| n = n \n —1, (34) 

for all values of n from — oo to + oo, and to agree with the factorial 
function when n is integral, that is, [1_ = 1, [2 = 1.2, |3 = 1.2.3, 
&c. We shall still call |n the factorial function, and (| n)~ l the inverse 
factorial. 


2 m 2 
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How, by the above 




(35) 


if x be positive, as we shall suppose throughout. We conclude that 

19 =r 

Further differentiations give 


V*° = p T , &o. (36) 

We, therefore, conclude that (| — l)~ l — 0, (|—2) -1 = o, &c., or that 
the inverse factorial function vanishes for all integral negative values 
of n. We therefore know the value of the inverse factorial for all 
integral values of the variable, and a rough curve can be readily 
drawn. Say 

V = (37) 


y being the ordinate, and n the abscissa. It has evidently a hump 
between n — 0 and n — X, is positive for all + valaes of n, asymptotic¬ 
ally tending to the n axis as n is increased, and is oscillatory on the 
other side of the origin. This is not demonstrative, but only highly 
probable so far. 


The Inverse Factorial Function. 

17. How seek an algebraical function with equidistant^ spaced 
roots on one side (either side) only of the origin. The function 




(38) 


vanishes at n = 1, 2, 3, &o., up to r. It has no other roots, and is 
positive when n is negative. Also, its value at n = 0 is 1. Similarly 
the function 


(i+^x+j) 




(39) 


vanishes at n = —1, —2, &c., up to —r; is unity at n — 0, and is 
positive when n is positive. These functions are identically the 
same as 

l-n+ —QQ— 2 ) + . + w( w —1)-(w-r+ 1) 

■ l! !« t ^ 

and 


n(n- 4-1) ( n (n -f-l) (n + 2) 

1 r n -j-—- 1 ---r * 


+ : 


> (w + l) . . . . ( n + r— 1) 
\ r 


(41), 
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where (40) corresponds to (38) and (41) to (39). At first sight, 
therefore, these functions might represent the inverse factorial, 
positive or negative, on making r infinite, for the value at the origin 
is correct, and the vanishing points are equidistantly spaced with 
unit step all the way to infinity on one side only of the origin. But 
something else happens when r is made infinite. The value of (38), 
by (40), becomes (l —l) % , meaning the binomial expansion in rising 
powers of the second 1. It is, therefore, zero for all positive and 
infinity for all negative values of n. Similarly, (39) becomes (l — l)~ w , 
which is zero for all negative and infinity for all positive values 
of n. That is, from vanishing at detached points, the functions vanish 
all the way between them as well. Besides, apart from this, we cannot 
have the value of j n correct when n is integral. 

We may, however, readily set the matter right. To get rid of the 
infinity on one side and vanishing all over on the other side of the 
origin, multiply the functions (38), (40) by r n and (39), (41) by r~' 1 . 
Take 

_i_ _ r -»(i_i)-« = r->‘( 1 + n )[ 1 +fj( y 1 + ' > ^) -> (42) 



We now satisfy all the requirements of the case, and when r is infinite 
make the inverse factorial curve (37) be a continuous carve from — co 
to 4-co, subject to (34), in agreement with the known values when n 
is positively integral, and harmonizing with the generalized dif¬ 
ferentiation in (33). 

Multiplying (42) and (43) together, we obtain 


1 1 

| n |—n 


(i-o 



sin nir 
mr 


( 41 ) 


The multiplication therefore brings all the equidistant roots into play, 
on both sides of the origin. 

This gives us the value of \—n in terms of | n. Only the values 
of |n from n = 0 to n = 1 need be calculated, since (34) or (44) gives 
all the rest. But if we take n — \ in (44) we obtain 


2 V I 2 ) — 2 


therefore |— (45) 

a fundamental result. We now know the value of when n is 

any integer, and this brings the matter down to the determination 
of \n from n — 0 to n — -jr, for which a formula may be used. 
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Interpretation of Vanishing Differential Coefficients. 
18. It should be noted that when we say that 


V 


x n 

n _ 

n 


= h 


(46) 


for all values of n? the right member is really x°/\o, and means 0 on 
the left, and 1 on the right side of the origin of x. That is, it is the 
limiting form of the function x n j\n, when n is infinitely small positive. 
It is convenient in the treatment of equations of the form F = Y f to 
have the function / zero up to a certain point, with consequently F 
also zero, and then begin to act. Similarly, the expression v^ 0 /l o > or 
yl or x~ 1 l\ — l, although it has .the value zero for all positive values 
of x, is infinite at the origin. But its total amount is finite, viz., 1. 
Imagine the unit amount of a quantity spread along an infinitely long 
line to become dll massed at the origin. Its linear density will, in 
the limit, be represented by, as (12) is derived from (11), 



cos mx dx. 


(47) 


It is zero except at x = 0. But its integral is still finite, being 
V°1 or 1. If we draw the curve y = x n j\ n, with n infinitely small, 
consisting of two straight lines, with a rounded corner, the curve 
derived from it by one differentiation will nearly represent the func¬ 
tion yl, being nearly all heaped up close to the origin, and of integral 
amount Similarly ,v 2 l means a double infinite point, v 3 l a triple 
infinite point, and so on. But it is the function v 1 that is most 
useful in connexion with differentiating operations, whilst the others 
are less prominent. 

But when n is taken to be infinitely small negative in y = x n j\n, 
then y drops from^oo to 1 near the origin, or the corner is turned the 
other way. That is, the function x n is unstable when n is zero. It 
is the difference of the curves y = x n with n infinitely small positive, 
and the same with n infinitely small negative, that makes the 
logarithmic ffunction\when infinitely magnified. But we should try 
to keep away from the logarithm in the algebraical treatment of 
operators. 


Connexion between the Factorial and Gamma Functions. 

18a. It will be seen by (42) that our factorial function is the gamma 
function of Euler somewhat modified and extended. Thus, when n is 
greater than — 1 we have 


jn = r(n + l), 


m 
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and this is also expressed by the definite integral 

p cr x x n _ ^ 

]o~T*: - 


(49) 


But when n is less than —1 we have the oscillatory curve of t\ 
inverse factorial, given by (44). We cannot use the definite integral 
to express |n when n is less than —1. The one-sided reckoning of 
the gamma function expressed in \n = F(w + l) is so exceedingly 
inconvenient m generalized differentiation that the factorial function 
had better be used constantly. For completeness and reference, we 
may add the general formula. Take the logarithm of (42) and 
arrange the terms suitably, and we obtain 



, i ~ n* „ nr ~ 

log I n= —nC +—S 2 ——S 3 +-- 

(50) 

where 

S “= 1+ ^.+^+-> 

(51) 

and 

C = Si — log r = 0*5772 ; 

(52) 


it being understood in (42) and (52) that r is made infinite. 

From (50) we may obtain a series for the inverse factorial in rising 
Powersoft. Thus, 

1 4,2 m3 

r = i+c»+(C 2 —s 2 )^-+(C 3 —3CS*+2S 3 £ 

\n A o 

+ (C 4 -6C 2 S 2 +8CS 3 +3S 2 2 -6S 4 )^+. (53) 

As before remarked, only the value of \n from n = 0 to \ needs to 
be calculated. Any number of special formulae for \n may be obtained 
from algebraical expansions involving this function. 


A Suggested Cosine Split. 

19. The above split of the function (sinmr)/^7r into l /|n and 
l/| —n, suggests other similar splits. In passing, one may be briefly 
noticed, the cosine split. Thus, take 


/(«) = (-!)"“ i+i i + 







(54) 


and let/( — n) be the same with the sign of n changed. Then, when 
r = co, we shall have 


= cos 7tx. 


(55) 
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By changing-^ to n +1 we find that 


f(n+ 1) = f(n) 


2n + 2r +1 
(2n + T){r—\) 


(56) 


so that when r = co, 

/(» + i) = ~7T ’ (50 

^ t 2 

or (»— 2 ) /(») = /(» — '!)• (58) 

When n = 0, (54) gives /(0) = 1. Then (58) gives f(n ) for any 

integral n. Thus /(l) = 2, /(2) = ,/(3) — -ft, /(4) = &o. 

And on the negative side we have /(-—■|) — 0, /( —1) = — 

/( —1|) = 0,/(— 2) = £,/(—21) — o,/(—3) — — 1 sS&c. The curve 
is similar to that of the inverse factorial, but with a much bigger 
hump on the positive side, near n = 1. But I have not found 
any use for this cosine split, and we may now return to the other 
one. 

The 'Exponential Theorem Generalized. 

20. Although we cannot, owing to its limited applicability, use 
Euler’s integral to express \n generally, we may employ it when 
found convenient, within its own range, and supplement the informa¬ 
tion it gives by other means. Thus, we know that 

f 00 e~ x x n 

l = -r— dx, (59) 

Jo \n 


when n is over —1. ISTow the indefinite integral may be exhibited in 
two different ways, say 



A>n +1 

jw + l 


+ 


C£ 7b 2 

1 _ I —LI_ L 

\n + 2 \n + 3 


= Wi, 


(60) 


in ascending powers of x multiplied by the exponential function, 
and by 




/yiTZ 1 /yt 72,-—~ 

+ j-- + |-- + • • 

\n — 1 \n — 2 


— w 2 . 


(61) 


These are true for all values of n. Subtracting (61) from (60) we 
see that the function w x -f- w 2 , or w say, must have the same value at 
any two finite limits we may choose for the integral. That is, the 
value of w is independent of the value of x. 

Or we may proceed thus, and determine the value. Let n be 
greater than — 1, and divide the integral (59) into two, one going 
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from 0 to a?, the other from x to oo. For the first use (60), since 
Wi = 0 when x = 0; anl for the second use (61), since w 2 = 0 when 
& = oo. We then get, by (59), 



1 = (WiH-^ 2 )j 


(62) 

or, which, is the same, 





/yiW 4" 2 /y>W“"l 

\ 1 9 P * • • • 4~ f- 4" I " 0 + • « 

In+2 \n — 1 \n — 2 


(63) 

~\n + \n + l + 



This is proved when n is greater than —1. But the change of n 
to n — l in the series on the right of (63) makes no alteration. We 
therefore conclude that the series expresses e x for all values of n. 

When n — 0 or any integer, positive or negative, we have the usual 
stopping series for e x . When n is fractional, we obtain semi-con¬ 
vergent series. Of course we obtain the whole series of forms by 
making n pass from 0 to 1. The most interesting case is that of 
n = J. This gives 


, 2x (2xf (2 ccY 

1+ ~3 + TT + ^J77 + ' 



-L-+±L. 

(2 a?) 2 (2 a?) 3 


1.3. 

J2P) 



(64) 


where the value of | we know to be ^| — \ = \'7f k , by (45), 

By means of this series we may pass from one to tlie other of the 
two forms of evaluation of Fresnel’s integrals, due to Knochenhauer 
and to Cauchy respectively, which are given in works on Physical 
Optics. 

21. The function called v.\ above we may obtain in a series of 
rising powers of a? without the exponential factor in the following 
manner:— 


t *—.r 

tVl = V ” 1 ——- — •— v~’(v + 

\n 
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To corroborate the method of getting (65) we may use (60). For 
this gives 

W 1 = 6“ A 'v'' {n+i) (l + V“ 1 + V“" 3 +_) = e-^v'“ (M+1) €' r = (v + l)~ ( ” +1) , (07) 


again, as in (65). How if we endeavour to express w 2 in a similar 
manner, we find that it will not work. But direct multiplication of 
the series in the brackets in (61) by the usual expansion of e~ x 
gives 


w 2 



jl—w+ 


n(n —l) 

“I* 


n(n—1 )(n—2) 

j3 






( n—1 )(n —2) 

+ JT 


(n—l)(n —2) (n — 3) f 



■ a a+1 f (w+l).ra (« + l).«(ro— l) 

' t "|w+l 1 ^ |2_ j3 

jjjm+s f («4 2)(» + l) (n+2)(n + l)n | 

+ 1 -ji- —[3 - + ’ ■ * * J 


or, which is the same, 

rM ntl-l q>1t + 1 

«?*= j —(1 —l)^ + j--(l-l)*” 1 *.... d“i“~r7(^ — l) n+1 "h* • • • 

\n ' |»—1 v |n4"l 


( 68 ) 


1 + 2 
FB+^ +1 >FEi 


(» + l)(» + 2) x n+3 

pr j^+3 + . 


(69) 


How the last line we know to express — w lt Therefore, by (62), 
we get 


rv,n q->n-T -1 

1= (i_i), + _ r _( 1 _ 1)B - i + .... + (i-i)^^...., (70) 


and this is the result we shall obtain by multiplying the series on the 
right of (63) by the usual expansion of e~ x . But (70) is only a 
special form of a more general formula that will appear later. We 
may use (44) to convert (70) to circular functions. 


A Bessel Function Generalized. 

22. The generalized expansion of e x may be at once applied to 
generalize other formulae. Thus, we know that the solution of 



\u — u 


(71) 
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in rising powers of x is 


T / \ , tv . w | *(/ 

lo («) = !+ 22+i^+^W+• • • ■ 


Put ar/4 = 2 /, and it becomes 


!«(*) = 1 + _ + _|_ + _|_ + . 


(72) 


(73) 


Let tbe differentiator djdy be called A. Then (73) may be at once 
written 


A" 1 A“ 2 A" 

Io(*) = l + ^+^- + |3- 


+ 


(74) 


There is nothing hypothetical about this. What is, however, 
entirely speculative in the absence of trial is whether it is legitimate 
to substitute the generalized exponential for the ordinary, and there¬ 
fore generalized for complete differentiations. But on trial it will be 
found to work. Thus, using (63) with A*"* 1 in place of x, we convert 


(74) to 


A A" (m+1) 

u = |-hi———b • 

jm I m +1 


a —(m—1) 

■+n~ T-+-- 


(75) 


writing u for the result; or, which means the same, 


}J m+l 


_(_£_L. 

(|m) 2 ^(|m+l) 2 ^ 


+ 


y w> 


(|m — l) 


2 +* 


(76) 


where, as before, y = oj 2 /4 in connexion with the Bessel function. 



if ^ = 2m. We have u m — Um^i for all values of m, and u 0 = I 0 ($). 

To test the validity when m is fractional, take m = | or n = 1, 
then, by (77), we obtain 


/y,3 /y,7 

- / %Aj %Aj tAj 

w = 


+-+-!+^f+—•?-+• 

X X 6 X X 1 


\ 


(78) 


This will by numerical calculation be found to give the same value as 
the series (72) for I 0 (a(). When x is as large as 10 the values agree 
to the fourth place by the convergent series in (78) alone, the semi- 
convergent series in (78) having a relatively small value. The value 
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I make to be 2815'71 by the ordinary formula, and 2815*75 by (78), 
including the semi-convergent part. But the last two figures are 
probably wrong, as there is a good deal of figuring involved in the 
calculation of both (72) and the convergent series in (78). 

When smaller values of x are taken, the numerical agreement per¬ 
sists as far as the initial convergency of the descending series permits, 
as in the case of the series (31), for example. Later on I will co¬ 
ordinate (78) and (72) with the descending formula (31). 

The companion formula to (78) is 



\x X d X 





^ 3 2 5 2 



(79) 


We might expect this to be a form of the oscillatory function J 0 (oj). 
But it is not. It represents the oscillatory companion to J 0 (a?), say 
G 0 (af), which maybe exhibited in an ascending* series of the whole 
powers of x 2 together with a logarithm, so standardized as to vanish 
at infinity. This function will appear later. The double series (79) 
occurs in Lord Bayleigh’s 4 Sound.’* The series (78) I have not come 
across. 

The Binomial Theorem Generalized. 

23. Let us next generalize the binomial theorem in a similar 
manner. We have 


(1 + aQ* _ 1 x x 3 

W \n |l \n — 1 |2 I n — 2 


(80) 


in ascending powers of x. Or 
(1 + aQ* = 1 . V" 1 V“ a 

| n | n | n — 1 | n — 2 


JL/Y! + Y!1_ + _ 

v* V n \ n —f 


(83) 


But in descending powers of x , which is the only other form gener¬ 
ally known, we have 


(a?+ l) M 


r~ + 




x Ql 


j n Jw — 1 (2 \n—2 






6 v 

V* 


(82) 


That we may use fhe generalized exponential we might infer from 
the two forms (80) and (82) being equivalent, combined with previous 


# See p. 154, vol. 2. That (79) is Gf 0 (#) is not explicitly shown there, but it may 
be readily deduced. 
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experience relating to the Bessel function. Using it in (82) we 
find 


(1 + aQ* _ |_V™ 


-+• 


which by immediate integration gives 


+ 1-” + • < 

\m — l 


(83) 


:+r 


, r n—m +1 


_ X n 

| n | m | n—m 1 \m —1 | n—m +1 


+ • 


+ 


[m + 1 \n — m —1 


+ (84) 

for all values of w. 

The case m = 0 or any integer is that of (82), and on ~ n is that 
of (80). The whole series of forms ranges between m = 0 and on = 1, 
because they recur. When on = ^ we have 




rpU + h, 


II K 


:i+r 


- n |« + i 


+ • 


+ T 


r + - 


(85) 


If in this we take n = 1, we find that the terms can be arranged 
in pairs, thus 


i-\-x 





05 f + OS 

Fill 


: + 


ad 4~ x f 

E 3 T 5 + 


(8G) 


or, which is the same, 


7r 

4 


(i+») = 


a* 1 + 


# + ® 1 

_T~-— 


x 2 -f x~ 2 t x ? ’ + x 3 
3.5 5.7 


(87) 


The best value of x is obviously 1. 

When two variables x and y are used in the binomial theorem, we 
have, using v for d/dx and A for d/dy, 


\n 


_ar* 

“ k + 


rpdft 1 

\oi — 1 


y+ 


rptl 2 ,i,2 

I —rF+ •• •• = 

I n — 2 [2 


A + .... 


yy-n 

1 + V/A 


1 + A/v 


— m A“ ?, e- rA 


( 88 ) 


We may use any of these forms. Selecting the last but one, and 
using the generalized exponential, we have 




f (?/v) M . Q/v)”- 1 . . (yv) 8,+1 , 

\ |m | m—1 |m + l 


•• } ; (89) 
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therefore 
(x + y) n __ y m x'‘ 


f,n—m ym—ly-n—m+I r^n—m—l 

+ 1— . r~-“TT + *•••+ T 


|n | m | n —m \m—1 | n — m-\-1 


|m +1 | n — m—1 

+ (90) 

where m may have any value, and x , y may be exchanged. 


Taylor's Theorem Generalized. 

24. We may also apply the generalized exponential to Taylor’s 
theorem for the expansion of a function in powers of the variable. 
For this theorem is expressed by 

/(» + *)==e*V(®)» (91) 

and, if this be true generally, irrespective of the wholeness of the 
differentiations, we must have 


r }> n 7,»+i 1 

/(*+'*) = {fT'-V +|S=i^ + -- } /W ' 


(92) 


Whether this is true for any function f(x), with the usnal limitations, 
I cannot say. There are probably other necessary limitations. 

As examples, take f(x) — 1. Then we obtain 


h n x~ n h n ~~ l x l ~‘ 
~ n]n— Y\T- 


\n 


-+ 


+■ 


h n+l x n 1 
n+1 I —n — 1 


+ 


(93) 


Here put c~hjx; then, by using (44), we have the result 
sin mr , sin(w-t-l)* 


(^ + 1)7 


_j sin (n — 1 )tt 
( n — 1)7 


where c is to be positive. When n = this reduces to 

i = ®+“*( a *+^)+* (« 5 +^)~ 


■, (94) 


(95) 


where a is written for c*. It is obviously right when a = 1. 

The formula (70) may be derived from (94) by the use of (44). 


Special Formulae for Factorials. 

25. The binomial generalization before given is, of course, a special 
case of (92), namely,/(a?) — x n j\n . It will be observed that the series 
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it gives may be convergent. Thus we may obtain convergent special 
formula! for in. Thus, take m = \, n — ■| in (84). We obtain 


( 1 +®) 


li 


(II ) 2 


i+ 


*+ or 1 3(a; 3 + x~'~) 3.7(x 3 + ai~ 3 ) 


+- 


--},( 96 ) 


5 5.9 * 5.9.13 

and when x = 1, we have the series 

(ff (I4)*=i+t(i-4(i-*(i-lf(I-.- - • (97) 

Similarly, m = £, n = £, as = 1, gives 

= 1—i(l—1(1 —tt( 1—if (!—xK 1 —-! (98) 


8r* 


81.4»([i) : 
and so on. 


Property of the Generalized Exponential. 

26. Notice that the operation performed upon the generalized 
e x reproduces it when m is integral, but gives an equivalent series 
when m is fractional. If, then, we take the special form of the 
ordinary stopping series for V w to work upon, we require to imagine 
that the zero terms are in their places, thus, 

e * = --'-+fr 2 +\ir 1 + 1 + x +^+ . ( 99 ) 

All terms before the 1 are zero, but not their rates of variation with 
x in the generalized sense, if we are to have harmony with the 
behaviour of the general form of e x . This is transcendental: and 
there is much that is transcendental in mathematics. 

The above generalizations are somewhat on one side of our subject 
of the treatment of operators, though suggested thereby. I propose 
to continue the main subject in a second paper. 




